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Abstract 

We consider reaction-diffusion equations of KPP type in one spatial di- 
mension, perturbed by a Fisher- Wright white noise, under the assumption 
of uniqueness in distribution. Examples include the randomly perturbed 
Fisher-KPP equations 

d t u = d 2 x u + u(l -u) + e^/u{l - ujW, (0.1) 

and 

d t u = d*u + u(l - u) + e^/uW, (0.2) 
where W = W(t, x) is a space-time white noise. 
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We prove the Brunet-Derrida conjecture that the speed of traveling 
fronts is asymptotically 

2-7r 2 |loge 2 r 2 (0.3) 
up to a factor of order (log | log e|)| log e|~ 3 . 

1 Randomly perturbed KPP and the Brunet- 
Derrida conjecture 

In this article we study randomly perturbed Kolmogorov-Petrovsky-Piscunov 
(KPP) equations, 

d t u = dlu + f(u) + ea(u)W £ > 0, xeR, u>0 (1.1) 

where W = W(t, x) is two-parameter white noise; / is assumed to be a Lipschitz 
function satisfying standard KPP conditions, 

/(0) = /(l) = 0; < f{u) < uf(0), u 6(0,1), (1.2) 

and in addition that for u > 1, 

f(u)<2-u, (1.3) 

We can and will rescale so that /'(0) = 1. We assume that tr 2 (u) is a Lipschitz 
function satisfying 

a 2 (u)<u (1.4) 
and for which there exist a* > and < u* < 1 such that 

<J 2 (u2) — u 2 (ui) > a*(u2 — ui), for < U\ < 112 < u* . (1-5) 

We will consider (jl.ip with initial data u(0,x) = uo{x) satisfying, for some 

x e R, 

uo(x) > 9 > 0, x < xq, and / uo(x)dx < 00 (1-6) 

J xo 

and contained in some subset C of the set C eX p of non-negative continuous func- 
tions / on R with f(x) < Ce' x ' for some C < 00, for which we know 

u {x)eC u(t,x) G C Vt > 0, P- a.s. (1.7) 

Weak uniqueness holds in C. (1-8) 

Key examples are 

f(u)=u(l- U ) (1.9) 

and 

a 2 (u) = u(l -u)l(u < 1) (1.10) 
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or 

(7 2 (u)=U (1.11) 

with uo(x) G [0,1] for (pTTP|) and u satisfying e~ x > u {x) for ([TTT]) . (fT§j) 
with p.lOp appears as the limit of the long range voter model and with (jl.lip 
as the limit of the long range contact process (see |MT95] V 

Note that (|1.3[) is not relevant for models such as (|1.10[) where < u(t, x) < 1 
for all time. But some condition on large u is needed in cases such as (jl.lip 
where fluctuations can take the solution above 1. 

We regard the stochastic partial differential equation (SPDE) (jl.ip as short- 
hand for the integral equation, 

u{t,x) = J G(0,y,t,x)u (y)dy + J J G(s,y,t,x)f(u(s,y))dsdy 
t 

G(s, y, t, x)a(u(s, y))W(ds, dy), (1.12) 

o 

where G(s, y, t, x) — G(t — s,y — x) is the heat kernel 

G(t, x) = (47rfT 1/2 cxp{-a; 2 /40, (1-13) 

and the white noise W is defined by specifying, for square integrable deter- 
ministic functions >p(s,y), that W(ip) := J J* °° ip(s,y)W(ds,dy) are a Gaussian 
family with mean zero and covariance 

E[W(ip)W(iP)}= J J <p( S ,y)ip(s,y)dsdy. (1.14) 

Here, and throughout, J fdx means the integral over the entire real line R. 
Solutions to (|1.12p are called mild solutions. See |Wal86] for the definition of 
the stochastic integral in (|1.12p . Readers unfamiliar with SPDE can think of 
the following system of ordinary stochastic differential equations on R; z , 

du i/n = [n 2 (u (i+1 y n -2u l/n +u (i _ 1 y n ) + f(u i / n )}dt+n 1/2 ea(u l/n )dB i / n , (1.15) 

where Bi/ n (t) are independent standard Brownian motions, and n is large. A 
corresponding evolution of functions on R is produced by connecting the points 
(i/n,Uif n ) and ((i + l)/n, it(i+i)/„) by straight lines, and (jl.ip is obtained in 
the weak limit as n — * oo. 

When e = 0, (Ti~Tjl with / of the form (fj~9|) is the standard KPP, or Fisher- 
KPP equation, introduced in 1937 by both Fisher ipj, and Komogorov, Petro- 
vskii, and Piscuinov KPP]. The basic facts in this case are: There is a one- 
parameter family F v of traveling front solutions F v (x — vt) with _F V decreasing, 
F v (x) — > 1 as x — > — oo, F v (x) — > as x — > oo and F v (x) ~ e~ lx for large x, 
with v = 7 + 7 _1 . For initial data Uq(x) = l(x < 0) we have convergence to the 
traveling front with minimal speed, 

v = 2 (1.16) 
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in the sense that 



lim sup \u(t,x) - F vo (x- m*(t))\ = (1.17) 

x£R 

where m*(t) defined by u(t, m*(t)) = 1/2 satisfies lim^oo m*(t)/t = vo- Further 
details about convergence of the KPP solution to the traveling front were given 
by McKean |McK75j and |McK76j , and Bramson [Bra78j and [Bra83j , among 
many others. 

When e > with initial data in C eX p satisfying (|1.6p . one has non- negative, 
continuous solutions, with a finite upper bound on the support 

r(t) = sup{x £ R : u(t, x) > 0}, (1.18) 

for t > 0. The process viewed from r(t), 

u(t,x) =u(t,x + r(t)) (1.19) 

should have a unique nondegenerate stationary solution. This is the random 
traveling front. One also expects t~ l r(t) to have a nonrandom limit, 

v e = lim r'r((). (1.20) 

t— *oo 

This was proved [MS95] in the case (|1.9[) . (jl.lOp , for sufficiently small e. They 
consider initial data < uq(x) < 1 such that uq(x) = 1 for x < t and uq = for 
x > r for some — oo<^<r<oo. Because in this case a (I) = 0, solutions stays 
within this class, and the result of |MS95j extends to any / and a satisfying in 
addition to our assumptions, that <r{u) = 1 + 0(y/l — u) as u f 1. 

We now make some comments to justify the form (|1.4j) . (|1.5p of the stochastic 
perturbation. The most important reason for taking 

<t(u) ~ y/u, < u « 1 (1.21) 

is that this is the type of correction seen when approximating the reaction- 
diffusion problems by microscopic particle models. For example, the reaction- 
diffusion equation with f(u) — u(l — it) was originally derived by Fisher as 
a model for the spread of an advantageous gene; the term u(l — u) in (|1.9p 
represents the frequency of mating between the individuals with and without the 
advantageous gene. If there is randomness in the mating, for example, if matings 
were successful with a certain probability, the variance of the random term is 
naturally proportional to u(l — u), and this leads to a term \J u(l — u)W(t, x). 

In this article we are primarily concerned with the asymptotics of v e as 
e — * in (|l.ip . It is not hard to see (for example, in (|1.9p . by taking expec- 
tations, and applying Jensen's inequality) that v c < vo. Recently, Brunet and 
Derrida |BD97j and jBDOlj (see also |KS98j . |PL99j ) have made the remarkable 
conjecture that as e — > 0, 



It is worth noting how enormous the correction is. For example, a naive Taylor 
expansion might suggest, since symmetry implies e = is a local maximum, 
that vo — v e ~ 0(e 2 ). The phenomenon was unexpected, and first observed 
through computer simulations of particle systems. It was not long before this 
was understood at the physical level as a consequence of the pulled nature of 
the fronts. Recall that in an evolution equation with traveling fronts between 
an unstable and a stable state, the front is said to be pulled if its asymptotic 
speed is the same as that of the linearization of the equation about the unstable 
state, and pushed if the speed is larger than that of the linearization (see [vSj ) . 
KPP equations have (marginally) pulled fronts. Because in pulled fronts the 
front speed is determined in the region where the density u is very small, in 
retrospect one should not be surprised that fluctuations there of order ^Ju, 
would have a dramatic effect on the front speed. Bramson |Bra78j also proved 
for the KPP equation with initial data uq(x) = l(x < 0), that, 

m*(*) = vo*+^log* + e>(l), (1.23) 

which is also supposed to be universal for pulled fronts [vSj . These behaviors all 
reflect the fact that in pulled fronts there is no spectral gap in the linearization 
around F Vo . 

The phenomenon (|1.22p has also been observed in systems where the variable 
u is forced to take discrete values, such as particle systems on the lattice with 
random walks and birth-death components. Here e 2 is the effective mass of a 
particle. 

Brunet and Derrida }BD97j conjectured that the front speeds in these sys- 
tems behave for small e like that of the solution of the cutoff KPP equation 

d t u = d 2 x u + u{l-u)\{u>e 2 ). (1.24) 



The idea is that when u < e 2 , u(l — u) < ey/u(l — u) and the noise term in 
(|1.1|) beats the creation term down to zero. Alternatively (ll.24|) can be thought 
of as a single particle cutoff. They then gave a nonrigorous argument that for 
small e, (|1.24p has travelling fronts with velocity 

7T 2 

Vcutoff — v - -j-j 575-. (1-25) 

|loge^ 

The argument for (|1.25[) . using matched asymptotics, is not difficult to make 
rigorous. It is known |DPK| . [BDLj that 

7T 2 1 

Vcutoff - v + j- ^ = 0(- 13). (1.26) 

|loge 2 | 2 |loge| d 

Implicit in our argument is a simpler proof (with a slightly worse correction of 

^/log_|loge|\\ 
|loge|3 ))• 

What was less clear was how to make rigorous the connection between either 
microscopic particle models or (|l.ip and (|1.24j) . Here we work with (|1.1|) as a 
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kind of canonical model system for the phenomenon (|1.22[) : In particular, the 
fact that particle models and (jl.ip both are expected to display this behaviour 
is perhaps the strongest motivation for the particular form of the noise (|1.5p . 

See [P] and references therein for a very comprehensive review of the physical 
aspects of the Brunet-Derrida theory. 

Conlon and Doering |CD04j recently obtained progress on (|1.22p by cou- 
pling (jl.9p . (jl.lOp to a contact process (see Liggett |Lig85| ), proving that for 
sufficiently small e, 

Clog | log e| 

v e > v r ,2—. (1.27) 

I log e| 2 

Very recently, [BDMM] have made a conjecture about the corrections to 
(jl.22p . Using a phenomenological argument, they propose 



- v + 



6tt 2 log I log e I 

— n„„.2i3 — ■ ^-28) 



loge 2 | 2 |loge 213 



The 6 on the right hand side is closely related to the 3 in (|1.23p . 

In this article we prove the Brunet-Derrida conjecture for models of the form 
(|1.1[) - (|1.3|) with a correction of the same order as conjectured in |BDMM| . 

However our understanding of the well-posedness of (jl.ip is not complete, 
and so a few comments are needed before we can state the result. 

Existence for is straightforward and can be obtained as the limit of its 
spatial discretization (| 1 . 1 5|) . Starting with non- negative initial data, one obtains 
in this way a non-negative solution, Holder a < 1/2 in space and j3 < 1/4 in 
time. Alternately, equations of the form (jl.ip can be obtained as appropriate 
limits of particle systems. Note that we are allowing solutions to have u > 1, 
which is slightly non-standard, in particular for models such as f(u) = cf 2 (u) = 
u{\ — u) where u is usually taken to be in [0, 1]. In terms of existence, this does 
not make any difference. 

On the other hand, uniqueness of (jl.ip is not known in our case because 
the coefficient in front of the noise in not Lipschitz. At the time of writing 
there is not even a consensus whether strong uniqueness should be true (for 
new results on strong uniqueness for stochastic partial differential equations 
with non-Lipschitz coefficients see |MPj . although they still do not cover the 
present case). Weak uniqueness means uniqueness of the martingale problem 
for (jl.ip with respect to the family of functionals f$(u) = exp{— J ucfidx}, <j> 
smooth, non-negative, with compact support, within the class of continuous non- 
negative solutions u, and in addition, the measurability of the Markov transition 
functions P S ^.)(A) = P(u(t,-) S A \ u(s,-) — u(-)). Weak uniqueness in 
particular implies the strong Markov property, which is one of our basic tools. 
Weak uniqueness can be obtained in special cases of (jl.ip using duality. An 
example is when a is of the form (jl.lpp or (jl.lip . The case (fL~9)) , (jl.lOp has 
an explicit dual particle system, described below, whose existence allows one 
in principle to compute the law for the stochastic partial differential equation. 
The case (TO]), (TrXTj) is self-dual (see [HT] V 
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What we really use about our solutions are the strong Markov property with 
respect to a family of hitting times, together with the comparison principle. 
Roughly the comparison principle for SPDE states that if u and v are solutions 
of 

d t u = dlu + f(x,u) + a(u)W, d t v = d 2 x v + g(x, v) + o-(v)W (1.29) 

on [0,T] with v(0,x) < m(0, x), and g(x,u) < f(x,u) for all x,u, then v(t,x) < 
u(i, x) for all t G [0, T] and x € R almost surely It is the analogue of mono- 
tonicity or attractiveness in particle systems. Examples of such theorems can 
be found in [Pa] . A simple variant of the above it that our v will, in addition, 
satisfy a Dirichlet condition v = on a set U such as x > vt or |a;| > vt + L. 
This is not a large leap, as one can think of it as the N — > oo limit of g = —N 
on U. So there is no surprise that the comparison continues to hold. There 
will be a few other twists and we will be a little more precise later. But the 
main point is that proofs of comparison theorems of this type require as input 
a strong uniqueness theorem. Hence they are not directly available to us. 

Now any solution we are really interested in will be the result of some ap- 
proximation scheme by systems, for example particle systems, for which the 
comparison principle is essentially obvious. Similarly, the strong Markov prop- 
erty will hold for such systems. And both are maintained under weak limits. 
So we could in principle just take a pragmatic approach and simply assume 
that our solution has the needed properties. Since this is a little cumbersome, 
instead we will state our results under the assumption of weak uniqueness. 

Note that weak uniqueness implies the existence of versions satisfying both 
the strong Markov property and comparison principle. That it implies the strong 
Markov property is well known. To obtain a version satisfying the comparison 
principle, construct a sequence of Lipschitz a^ n \u) converging uniformly to 
o~{u). The corresponding equations have strong uniqueness and therefore the 
comparison principle. It is not hard to check that such sequences are tight and 
it is easy to see that the comparison principle continues to hold in the weak 
limits. Note that all our results are statements in distribution. It is therefore 
always enough to work with appropriate versions of our process, and therefore 
weak uniqueness is sufficient. 

We can now state the main theorem. Let u(t, x) be the solution of (jl.ip and 
r(t) be as in (|1.18p . For initial data in C exp satisfying (|1.6[) let 

v e = limsupt _1 r(i). (1.30) 

t — >oo 

For initial data satisfying u(0, x) > 9 > 0, x < 0, u(0, x) — 0, x > 0, let 

v e = ]imin£r l r(t). (1.31) 

t — >oo 

Let a(a) be the largest a such that 

(1 - a)ul(u < a) < f(u). (1.32) 

Note that a{a) > if a > from the assumptions on /. For example, if 
/"(0) > — oo we have a(a) = Ca for some C < oo. 
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Theorem 1.1. Assume Uq(x) £ C satisfying il.6}) , |J.7[ ), U.8\) . Then there 
exists £q > such that for all e < eg, 



X e > v 
v e < v 



27r 2 [91og|loge| -loga(|loge| 



log e 2 1 2 



loge 2 |3 



loge 2 



87r 2 log | loge | 
I loge 2 | 3 



(1.33) 
(1.34) 



Remark on cq. Since the phenomena is observed in particle simulations, it 
is worthwhile to ask whether the mathematical result can be proven with eo 
of a size approachable by computation. In fact, computations with effectively 
N = 10 10 particles are commonplace at the time of writing. In case fj 1 . 9[) . 
(|1.10|) . we can check that our method works with eo = e~ . Since N particles 
corresponds to a correction of size e = TV -1 / 2 the mathematical result covers 
the typical regime of computations. 

Remark on duality for the case of il.l]) with f(u) = o~{u) = u(l — u). Let 
Xi(t) : i = 1, . . . , N(t) < oo be a system of independent Brownian motions with 
generators <9 2 . Each particle splits in two at rate 1, and pairs of particles coalesce 
at exponential rate e 2 during their intersection local time. The generator is 

a/( X ) = A </« + E(/( x +) - /(*)) + e2 E **<-i (/( x ^) - ( L35 ) 

i i i>j 

where is the configuration obtained from x by replacing Xi by two particles 
at the same location, and x~ is the configuration obtained from x by removing 
Xi. We have the duality relation |Shi88j . 



E 



H(i- U (t, Xi (o))) 



= E 



11(1 - u(0, 



(1.36) 



where the expectation is taken over independent u and Xi. Among other things, 
(jl.36p gives us an expression for the moments of it, providing the weak unique- 
ness. 

One can also deduce from the result about the wavespeed in the random KPP 
results about the wavespeed in the dual process. Suppose we start our branching 
and coalescing system with one particle at 0, and let L(t), R(t) denote the 
positions of the leftmost and rightmost particles in the system at time t. Take 
uq(x) = l(x < 0). The duality relation, together with the natural reflection 
symmetry and spatial homogeneity, give P(L(t) < —x) = P(R(t) > x) = 
E[u(t,x)}, and Theorem 1 1 . 1 1 then translates to 



lim = - lim £ -1 L(t) 



(1.37) 



Here then is another example of the Brunet-Derrida theory: The branching- 
coalescing Brownian motions model possesses two invariant measures. The sta- 
ble one is a Poisson point process with intensity e~ 2 , and the unstable one 
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consists of no particles at all. On a large scale we see the first invading the 
second at linear speed v £ . If we introduce a phase variable in [0, 1] so that 
corresponds to the unstable phase and 1 corresponds to the stable phase, then 
the effective particle mass is e 2 , as predicted. 

Finally we comment on the structure of the paper. To make the arguments 
leading to Theorem 11.11 more transparent, in the next section we sketch the logic 
of the proof, assuming the main technical lemmas, which are then left for later 
sections, and assuming as well that the necessary manipulations of the SPDEs 
can be performed. We then prove the validity of these manipulations in Section 

El 



2 Outline of the proof 

2.1 Comparison equation 

The general idea behind our proof of Theorem 11.11 is to compare the stochastic 
KPP evolution (TEl]) to: 

J d t g = d 2 x g + f(g), x<vt, 

\ g(t,x) =0, x > vt. [ ' 

We search for the v — v com for which there exists a traveling front solution 

q{x, t) = F(x - v com t) (2.2) 

with lmxr^-oo F(x) = 1 and 

F'(0) = -e 2 . (2.3) 

The problem (|2.1[) - (|2.3p is our replacement for Brunet and Derrida's com- 
parison equation (|1.24[) . The idea is that the solution will have a mass of C(e 2 ) 
within a distance 0(1) of x = v com t. Here e is a small parameter which docs 
not necessarily have to be related to the perturbation parameter e in (jl.ip . In 
fact, it will be convenient to take e to be slightly larger or smaller than e. But 
if e = e, the mass 0(e 2 ) is the critical mass which can survive in the stochastic 
equation when u is small. Heuristically, this will provide a consistent strategy 
for a stochastic traveling front in (II. ip to propagate. 

To determine the resulting v = v com (e 2 ), let x(i) = F(—t) and note that 
the problem is equivalent to that of finding the v such that the solution of the 
ordinary differential equation 

x" = vx' - /(x), x(0) = 0, x'(0) = e 2 (2.4) 

has x(oo) = 1. Let y = x'. In the phase plane of 

x' = y (2.5) 
y' = vy - /(x) 
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there is an unstable node at (x, y) = (0, 0) and a saddle point at (x, y) = (1, 0), 
joined by a separatrix solution (x(t),y(t)), — oo < t < oo, with x(— oo) = 0, 
x(oo) = 1. For v > 2, the linearization xj in = yi; n , y[ in = vyi in — xi;„ about (0, 0) 
has distinct positive eigenvalues which merge as v J. 2, then split into a complex 
pair for v < 2. For v > 2, the separatrix corresponds to an exponentially 
decaying traveling front solutions of the (nonrandom) KPP equation. For v < 2, 
(0, 0) is a spiral source, the x > corresponding in the same way to a traveling 
front solutions of (|2.ip . The separatrix enters the region (x, y) e (0, 1) x (0, oo) 
at (x, y) = (0,e 2 (v)), and problem (|2.ip - (|2.3p is now seen to be equivalent to 
computing the inverse function v(e 2 ). 

Proposition 2.1. Let v com = v com (e 2 ) be the solution of l£.i|) - fO)) . There 
exist £o > such that for \e\ < So, 

TT 2 TT 2 

(| loge 2 | - 3 log | log£ 2 | - Iogo(| loge 2 |- 3 ) - 2) 2 " Vcom " (| l og£ 2| + 3) 2 ' 

(2-6) 

If a(a) = a then e can be taken to be e 8 . 

First we give the heuristic idea of the proof. It is not hard to see that v(e 2 ) 
is monotone decreasing in e with v J. 2 as e 2 J. 0. The linearization about (0, 0) 
has explicit solution 

xjia(t) = eH- 1 ' 2 exp{(l - |) t} sin(S 1/2 t) (2.7) 

where 

S = 2- v. (2.8) 

and 5 = 6—^. Let Q be the smallest t > with 

xii„(9) = 1. (2.9) 

Note that ~ | log e 2 1 . Pretending the linearization is meaningful globally, one 
would want 

xn„(©) = 0. (2.10) 
(EU) and (i2TTU|) become, with d = 1, C 2 = (1 - |), 

e 2 <T 1/2 exp{(l - |) 9} sin(^/ 2 e) = d (2.11) 
e 2 exp{(l- |)e}(-cos(5 1 / 2 6)) = C 2 . 

(|2.1ip gives a nonlinear equation for 8 in terms of e from which it is simple to 
obtain estimates like (I2.6p . The only difference in the rigorous proof is that we 
will use sub- and super-solutions to get (|2.1ip . but with slightly worse C\ and 
C 2 . 

Proof of Proposition \2.1\ First of all note that v = v com depends monotonically 
on /: If fi(u) < / 2 (w) for all u then the corresponding v(/i) < v(/ 2 ). 
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Upper bound. Consider (|2.ip - (|2.3p with / replaced by 

/(u) = {; i ! uG[ i ) ' 1] (2.i2) 

The corresponding v = v com is larger as / > /. Call 5 = 2 — v and assume 
momentarily that 5 < 0.2. The solution to 

x"=vx'-/(x), (2.13) 

can be computed explicitly. For t < = inf{£ > : x(0) = 1}, x (£) = xi;„(£) 
from (|2.7p . Now we consider the phase plane (x, x'). One checks that the 
linearization around the saddle (x, x') = (2,0) has stable direction (w,l) and 
unstable direction (A, 1) where, 

w{8) = 1 - | - (2 - 5) 1 ' 2 , A =l-| + (2- 5) 1/2 . (2.14) 

A separatrix solution (x, x') joins (0,e 2 ) to (2,0). It must coincide with the 
stable line y = vox — 2m in the region x > 1, because the equation is linear 
there. So in order for (x(0),x'(0)) to lie on the separatrix we must have 

4 n (0) = wx lin (0) - 2w, (2.15) 

which is equivalent to (|2.1ip . with C\ = 1 and C 2 = — zu. 

Lower bound. Let a = a(\ loge 2 | -3 ) from ()1.32p and a = 1 — 2| loge 2 | -3 and 
define 

7, s / aw if m < a/2 

/(m) = \ a(f-tt) ifu>a/2. (2 - 16) 

Consider the problem l|2.ip - (|2.3p with / replaced by /. The corresponding 
v = v(/) is smaller than v com because / < /, Call 5 = 2 — v. From the upper 
bound we know that for sufficiently small e , S = 5 — ^- > 0, in which case the 
solution of 

x" - vx' + /(x) = (2.17) 

with x(0) = and x'(0) = e 2 is x(t) = x Un (<) from for £ < when 

xnn(0) = ot/2. In order to lie on the separatrix joining the unstable fixed point 
(x, x') = (0, e 2 ) to the saddle point (a, 0), we must have (x(0),x'(0)) on the 
stable line y = (1 — y2 — a)(x — a). This gives (|2.1ip with C\ — a/2 and 
C 2 = (v/a - ^- 1 - |)a/2. 

Proof of i2.6]) . Assume e n is sufficiently small that S < 1 and. Dividing the 
two equations in (|2~TT1) gives - tan(<5 1 / 2 0) = <5 1 / 2 (Ci/C 2 ). Let = S-^ 2 n~f3. 
Since < tana; < \/2x if x < tt/4 the solution has < C l /(V2C 2 ) < (3 < 

7r/4. Now the second equation of (|2.1ip gives, 

6 (1 - I)" 2 = tt 2 (| loge 2 | + logC 2 + I logcos^ 1 / 2 /?)! + p)~ 2 (2.18) 

To get a lower bound on v com = 2—5, drop the non-negative terms | logcos(<5 1 / 2 /3)| 
and (3 from the right hand side and note that 5 < 5 (1 — f) 2 and logC 2 = 
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\og{(V2~~ti - l)a/2} > -|loga| - 3 log | loge 2 | - 21og2. To get an upper 
bound, note first that if 6q is sufficiently small, then from the lower bound we 
have just described, 5 < 1. Since (3 < 7r/4, we then have | logcos((5 1 / 2 /3)| < 1. 

Also logC2 = log (V2-6- 1 - 6/2) < 0. Finally, if we take £0 sufficiently 
small, then the lower bound we just proved gives S < 10| log e 2 1 ~ 2 so that 
5(1 - !)~ 2 > S(l + 10| loge 2 !" 2 )- 1 and then ((2~T8)l gives 

6 > tt 2 (| loge 2 | + 2)- 2 - 100| logs 2 !" 4 (2.19) 

which gives (|2.6p for £0 small enough. □ 

2.2 Upper bound 

Consider u satisfying (jl.ip with initial data 

u(0,x) = F(x) (2.20) 

where g(t, x) = F(x — vt) is a traveling front solution of 

/ d t g = dlg + f(g), x<vt, . . 

\ g(t,x)=0, x>vt, { > 

with / as in (f2~T2|) . 

v = v com (£ 2 ) + |loge|~ 3 (2.22) 

and 

- F'(0) = £ 2 = 7e 2 (2.23) 

with < 7 « 1 to be chosen. 

F is a modified version of a traveling front from the comparison problem 
(|2.1[) - (|2.3[) . with a larger / > /, a slightly larger speed, and lying slightly above 
the separatrix connecting (0,0) to (2,0). There is some convenience in using 
/ instead of /. It is convex. Also, some things are explicitly computable. For 
example, 

( 0, x>0; 

F(x) = < xiin(-x), -6 < x < 0; (2.24) 
[ K\e- x( - e+x ^ +2- (l + /t)e- CT ( 0+:E ), a;<-e, 

where k — k(e 2 ) is chosen such that F and F' are continuous at x — — and 
xun and 6 are defined in (|2.7|) - (|2.10[) . Keep in mind that these all depend on 
e though the dependence is not written explicitly. One can check that k ~ 
(1 — 2 _1 / 2 )| log£ 2 |~ 2 . Note also that the modification of the speed in (|2.22|) is 
smaller than the ) error terms in the main result, Theorem ll.il 

Fix a positive integer T and an L > and consider the hitting time 

£ = inf {t e [0, T] : u(t, x) > F(x - vt - L) for some x £ R} (2.25) 
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We run u up to time £, and then restart with new initial data 

F(x - v£ - L), (2.26) 

a shift of L from the original comparison front. By the strong Markov property 
and the comparison theorem (Proposition 13. ip , we obtain an upper bound u 
on the original solution of with initial data F. Repeating the process, we 
inductively define a sequence of stopping times £,•_(- 1 € £j + T], and an upper 
bound G for all time on the solution u with initial data F(x). u satisfies (|l.ip 
on (^,^i_|_i) with initial data u(£i,x) = F{x — v£, — Li). 
Suppose we can show that 

P(£<T)< 1/2. (2.27) 

By the law of large numbers, the speedup of the front of u over that of u is by a 
factor L/E[g\. But from (|2~2Tj) . E[£] > T/2. We obtain in this way, using 
(|2.22p an upper bound on v e defined in (jl.30p . 

v e < v + 2T~ 1 L 

< 2-7r 2 (loge- 2 + 3) _2 + |loge|- 3 + 2T- 1 L. (2.28) 

If we choose 

T^L^Iloger 3 (2.29) 
we obtain the upper bound (|1.34p for initial data bounded above by F. 

There is a tradeoff between T and L. Large L in principle makes (|2.27p 
easier, because F(x — vt — L) is increasing in L. But then (|2.29p forces us to 
choose T large, and it becomes harder to control u on the long time interval to 
obtain ([2~27]l . 

For more general initial data, satisfying only J Q uo{x)dx < oo, we can use 
the fact that at any time t > 0, r(t) < oo a.s. This is proved in the special 
case <J 2 (u) = u in [MP92] . but it is well-known that the method can be adapted 
without too much work to cover the present situation. We do not give details 
here. This means that we can bound u(t, ■) by a shift of F, and obtain the upper 
bound (|1.34p as before. 

This reduces the upper bound to (|2.27p . The main idea to prove (|2.27p is to 
split the solution u(t, x) of (II. ip with initial data F into 

u = v + w (2.30) 

where v(t, x) is the mass which does not cross x = vt; 

d t v = dlv + f(v) + ea(v)Wi, x<vt, , , 

v(t,x) = 0, x > vt, [ 1 

with v(0, x) = F(x), and w > is the rest. Wi will be another space-time white 
noise. As usual, the SPDE is interpreted in the mild sense; 

v(t,x) = [ G v (0,y,t,x)F(y)dy + J J G v (s,y,t,x)f(v(s,y))dsdy 

G v (s, y, t, x)a(v(s, y))Wi (dsdy) (2.32) 
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where G v (s,y,t,x) is the sub-probability density for a Brownian motion with 
generator d 2 , starting at y at time s, to end at x at time t > s never having 
entered the region {z > vit} for times s < u <t. In Proposition ^. II in Section[3] 
it is shown that we can find a probability space on which such a splitting holds. 

We expect the solution v(t, x) of (|2.31|) to remain close to the solution g(t, x) 
of the deterministic comparison equation (|2.ip with the same initial data F(x). 
Because it is a subsolution of (|2.2ip we have g(t, x) < F(x — vt). If L is large 
enough, we can therefore expect v not to hit F(x — vt — L) for some time. 

The key point now is that if F'(0) « e 2 then w is so negligible that u — v+w 
does not hit F(x — vt — L) for some time either. To prove this, we need a better 
way to represent w. One can also view the Dirichlct boundary condition in 
(|2.3ip as a removal of mass. Let A(t) be the mass which is removed at the 
boundary x — vs in (|2.31[) during the time interval < s < t. Then we have 
another representation for v satisfying (|2.3ip or (|2.32p : 

d t v = dlv + f{v) + ea(v)Wj. - 8 x - vt A. (2.33) 

We would like to write an equation for w, with a new white noise W2, 
independent of W\ . If W\ and W2 are independent white noises, then 

a t Wi + a 2 W 2 = \J<jI + <jI W (2.34) 

where W is a white noise. Hence the equation for w should read, 

d t w = d 2 x w + f(v + w)- f{v) + eaW 2 + S x - vt A, (2.35) 

with initial data w(0,x) = 0, where 

a(t, x, w ) = ^a 2 (v(t,x)+w)-a 2 (v(t,x)). (2.36) 

But this is only reasonable as long as a 2 (v(t,x) + w) — a 2 (v(t,x)) remains 
non-negative. In Proposition 13.11 of Section O it is shown that there exists a 
probability space on which there are white noises W\ and W 2 for which (|2.35p 
holds, up to a stopping time 

r = inf{f > : 3x, a(t, x, w(t, x)) = 0, w(t, x) ^ 0}, (2.37) 

after which the desired noise coefficient a might cease to make sense. 

By the comparison theorem, and since f(v + w) — f(v) < \\f\\-Li p w, up to 
time r we have u — v — w < w almost surely, where 

d t w = d 2 w + WfWupW + ea(w)W 2 + 6 {x=vt} A. (2.38) 

As long as 

a 2 (w) > a*w, (2.39) 

this is basically a superprocess with an injection of mass at {x — vt}. The 
critical input of mass in such an equation is easily calculated to be 0{e 2 ). In 
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other words, if the rate of mass entering is o(e 2 ), then it is being killed by the 
noise in time 0(1) with very high probability. And it suffices to show just that 
the expected incoming mass E[A(t + 1) — A(t)\ is o(e 2 ). 

To get such a bound, note that by comparison v < v, the solution of 

d t v = dlv + f{v) + ea(v)W, x<vt (2.40) 

with v = for x > vt. Take expectation in (12.40)) and use the concavity of / to 
see that E[v] is a subsolution of (|2.2ip . In particular, 

E[v(t,x)]<F(x-vt). (2.41) 

This can be translated into a bound on the expected rate of incoming mass A(t) 
as follows. Taking expectation in (|2.33|) . 

E[A(t + l)-A(t)] = fq(t,y,t,t+l)E[v t (y)}dy (2.42) 

t+i 

q(s, y,s,t+ l)E[f(v(s, y))]dsdy 

where q(s,y,u,t) = P StV {3r G (it, t] : B r > vr) for a Brownian motion B r with 
g_enerator d 2 . Using E[f(v)} < E[f(v)] < f(E[v]) < f(p), and that v Q (y) = 
F(y), we see that 

E[A(t+l)-A(t)} < J q(t,y,t,t+l)F(y-vt)dy (2.43) 



t+i 



q(s, y,s,t+ l)f(F(y - vs))dsdy. 

Now F(x — vt) is a traveling front solution of (|2.ip with / instead of /. The 
rate of mass removal at the boundary is proportional to the slope e 2 at the 
boundary, and hence there is a C j2.44t < 00 sucn that, 

E[A(t + l)-A(t)] <C M £ 2 . (2.44) 

The only difficulty is maintaining (|2.39p . By (|1.5|) it holds as long as w + iD < 
u*. What we will do is obtain an a priori estimate that v < u* /2 in a strip 
vt — M < x < vt. This is reasonable since we know that i? is close to p, which 
is of 0{e 2 ) there. If 

M = M{u*) (2.45) 

is chosen sufficiently large, we can then iteratively show that w < u* /2, and 
furthermore that it does not support the complement of a strip vt—1 < x < vt+1 
around our proposed front. This then provides us with sufficient noise to show 
that w is negligible there as well. 

To fix 7 and T, let us explain very briefly the iterative procedure. Take T 
to be an integer and divide up the time interval [0, T] into intervals of length 
1. The mass arriving in [n, n + 1) and evolving according to (|2.38p . is bounded 
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by (|2.44j) . It dies before time n + 2 with probability at least 1 — C07 where 
c o = co(a*, ||/||Lip)' So this happens for every n = 0, 1, 2, . . . , T — 1, with 
probability at least 

1 - c Q jT. (2.46) 
In order to have the probability in (|2.46[) greater than |, we thus take 

7 = \^ X T- 1 . (2.47) 
To fix all our constants we note that if 

L = |loge|+log|loge| (2.48) 

then we have 

F(x -L) > F(x) + 3\e- Xx . (2.49) 

By (|2T2"9"j) we need 

r=L|loge| 4 J. (2.50) 

We have explained how the upper bound is a consequence of the following lem- 
mas. 

Lemma 2.2. Let v be the solution of \2.31)) with initial data {2.2$ . Let g 
be the solution of i2.21\) with the same initial data. Let j,L,T,M be as in 
(24^-(EM>- Then 

P(3te [0,T] :v(t,x) > g(t,x) + 3\e- x( - x - vt) for someiGKj < 1/16. 

(2.51) 

Lemma 2.3. Under the same conditions as in Lemma \2.2\ 

P(3t e [0,T] : v(t,x) > u*/2 for some x G (vt-M,vt)) < 1/8. (2.52) 

Lemma 2.4. Under the same conditions as in Lemma \2.2l let w be the solution 
of JUMP where A is defined in 

P(3te[0,T}: sup w(t,x)>0 or sup w(t, x) > u* /2) < 1/4. 

a:e[vt-l,vt+l] c x£[vt-l,vt+l] 

(2.53) 

2.3 Lower bound 

The proof of the lower bound uses a more standard method; comparison to 
oriented percolation [BraDurr88j. Similar arguments were used by Conlon and 
Doering |CD04j to prove their lower bound. The improvement here comes from 
the use of the special comparison front from (|2.17[) and refined large deviation 
estimates. 

For v, L > let G v ^L(s,y,t,x) be the sub-probability density in x at time 
t for a Brownian motion B Ul s < u < t with generator starting at B s = y, 
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killed if it enters the region \z\ > vu + L, s < u < t. If g(t, x) is a given function, 
let 



t 

2 



(Qg)(t,x) = J J G vL (s,y,t,x)g(s,y)dsdy (2.54) 

Note that this makes sense since we are in one dimension. The lower bound is 
based on the following simple lemma about the deterministic equation. Let 

/(„) = / £-\]° ge \~l\" , l U - a/ ^ (2.55) 
- v ' \ (1 — | log e| )(§-■") if u > a/2. v ; 

From the definition (ll.32p of a, we have 

/(«) < /(«) whenever f(u) > 0. (2.56) 

Let 

e 2 = 7 e 2 (2.57) 
with 1 « 7 and v = v com (e) as in (12. ip - (|2.3[) . 

Lemma 2.5. TTiere exist €q > 0, C jg < oo, < L < | log e | , and f? Q (a;) 

supported on [—L,L] with < £> (x) < a(|loge 2 | -3 ) smc/i £/iai if e < eo, 7 > 
C \2.5(h I l°g e| 10 and r > 1, t/ie solution g(t,x) on < t < 1 0/ 

<9 t £ = <9 2 £ + /(g), |z| < L + vt , , 

g(t,x)=0, \x\>L + vt, 1 M) 

with £>(0, x) = g Q (x), and f as in \2.55\) . satisfies, for v' = v + | loge|~ 3 , 



g(l,x)-reyjgg(l,x) > Q (x-v') x G [v' - L, V + L]. (2.59) 

Proof. We follow the notation and construction from the proof of the lower 
bound of Proposition ^. II Let a, a, x and O all be as in the proof of the lower 
bound of Proposition 12.11 We claim that (|2.59[) holds with L = and 

To prove this, note that the solution g(t, x) of (|2 . 58[) satisfies git, x) < q(t, x) — 
e *l lo s £ l 3 min{x(L + vt~ x), x(L)} since gJx) < x(L — x), x € R, and q(t, x) is 
a supersolution of (|2.58|) . and g[t,x) > q(t,x) where 

f e*l loge r 3 x(L + vt- \x\) vt<\x\<L + vt 
q(t,x) = i e 4 l lo « e l" 3 x(L) |x| < vt (2.61) 

[0 jxj > vt + L 

since g(f , ar) is a subsolution of (I2.58| with the same initial data. So 

£(1, x) - g Q {x - v') > q(l, x) - £ Q (x - v') (2.62) 



> 



x lin (i + v- |a;|) -x Ua (i + v' - |a;|) |z| e (v,v' + i], 

(e iio ge r 3 _ 1)xi . n(L) N < Vj 
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Here we used that x(t) = xii„(t) for t < L. For |x| € (v, v' + L] we now use that 
if 8 < 1/2 then x| in > £ 2 e5( i+v -M) there and v - v' = | loge|~ 3 to get a lower 
bound. For \x\ > v we use e I lo s e I ~ 3 — 1 > | log e|~ 3 and x lin (i) = a/2. This 
gives 

£ (l |a 0-4,(*-v)>| | lo ^j_ 3a/2 ' |a; |< l v ; (2-63) 

From the explicit form of xii n , and using G v .l < G, as long as exp{| loge|~ 3 } < 2 
and 8 < 1, 

1 vs — ' M_x ' 2 

< 2 / / , 6 7t~~~~~t min{ £ 2 r 1 / 2 e vs -^+ L , a /2}d S ^ 

Jo J-oo 4tt(1 - s) 

< min{4£ 2 5- 1 / 2 e-( x - E '- v ),a}. (2.64) 
Hence one can check that g(l, x) — g Q (a; — v') > re^/Q^ as long as 



r < min{67 1/2 |loge|- 7/2 ,e- 1 |loge|-Va(|loge|" 3 )/2}- (2.65) 

Since Qq>Qg we are done. □ 

From the comparison theorem (Proposition 13 . lj ) . we can construct a proba- 
bility space on which the solution of 

d t u = dlu + l(u) + ea(u)W, \x\<L + vt , , 

tt(t,aj)=0, |x|>L + vt. 1 ' 

gives an almost sure lower bound for the solution of (|1.1[) with the same initial 
data. 

Suppose we start (|2.66p with g (x) from Lemma [2.51 The idea is that the 
solution u will stay close to g up to time 1. To see how close, let us make a very 
rough argument. Since / has Lipschitz constant 1 one expects that for times of 
\u(ti x ) ~ Q(t> x )\ is controlled by something like e\ZL, v (t,x)\ where 

Z L ,vfo*) = J J^e t - s G v . L {s 1 y,t 1 x)G{u{ S ,y))W{dsdy). (2.67) 

The actual bound is somewhat more complicated, but it amounts to the same 
thing. Recall that <J 2 {u) < u. By concavity of /, E[u] < g. So 

E[Zl v (l,x)]<eGg(l,x). (2.68) 

Things are tight in the region close to the front {x = vt + L} where g = 0{e 2 ). 
Here Qg(l,x) = 0(e 2 ) as well. Hence the fluctuations of u — g are of order 
ee = je 2 there. This is the reasoning behing (|2.59[) . The key point of the 
following refined large deviation estimate is that it shows that the fluctuations 
near the front are of 0(je 2 ) instead of the 0(e) one would obtain naively. 
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Lemma 2.6. Let g, L be as in Lemma \2.5[ There exists C \2.(j{fa < 00 such 
that for all < r < e , 

p(Bx:u(l,x) <Q(l,x)-re^QQ{l,x)) < 4£exp{-C^r 2 }. (2.69) 

Note the factor L on the right hand side is because the large deviations are 
done on space intervals of size 1, and then summed over the width of g(l, x). If 
we take 

r 2 =C B ia[llogel 3 + 21oglloge]+log4] 1 (2.70) 

the right hand side is less than | loge| -1 e~l loge l 3 . We conclude that if we start 
(|2~66)) with g Q (x), then 

P(u(l,x) > g Q (x-v), x G R) > 1 - I loge|- 1 e~ |lose|3 . (2.71) 
Now we can ask for this to happen T — | loge| times, to obtain 

P(u(T,x) > g (x~vT), ieR)>l-e-l losf l S =j,. (2.72) 

By symmetry we also have 

P(u(T, x) > g Q {x + vT), x G R) > p. (2.73) 

This allows us to compare the system to a 2-dependent oriented percolation. 
Let C = {(to, n) € Z 2 : to + n is even, m > 0}. Let C denote the set of directed 
bonds (to, n) — ► (to + 1, n + 1) or (m, n) — * (to + l,n — 1). Let X&, 6 G £ be 
random variables taking values in {0, 1}. Assume that for all b 6 £, 

P(X h = 1) > p. (2.74) 

Also assume that Xb and AT^' are independent if the lattice distance between b 
to b' is strictly larger than 2. If toi < 7712 then we say (toi,71i) — > (7712,^2) if 
they are joined by a sequence of directed bonds bi € £ with = 1. Let iS 
denote the subset of (to, n) G C such that (0, x) — > (to, n) for some x < 0. 

Lemma 2.7. Suppose that 

p>l~e~ 10 . (2.75) 
Then, with probability 1, /or 6u£ finitely many to, 

A m d = max{n : (m,n) G 5} > (1 - 10| log(l - p)| _1 )m (2.76) 

Proof. Using the standard contour counting argument one obtains 

P(N m < to(1 - 6)) < E n > mS 4 2 " +m (l - P) n/2 - (2-77) 

(see |CD04j . Lemma 3.6, for a complete proof of (|2.77p . The only difference is 
the exponent n/2 on the (1 — p) which comes from the 2-dependence.) Take 
8 = 1 1 log(l — If we assume (|2 . 75[) then it is not hard to check that the 

right hand side is less than 2e" M and the result follows from Borel-Cantelli. □ 
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Remark. (|2.76[) is very far from optimal: It is known (see [D], |GPj ) that for 
p close to 1, the speed of oriented percolation, linim^^ ^ = 1 — 0(1 — p). If 
one uses the stronger result, one can check for the case (|1.9[) . (|1.10|) . the main 
result holds with e = e 11 • 

Proof of the lower bound. Start with initial data uo in C exp satisfying 
(|1.6p . Without loss of generality, xq = 0. As long as 9 > a we have Uo(x) > 

Z) n =-oo£o( X ~ W<3 Say ^(m,n)^(m+l,n±l) = 1 if "K.^ nT),u({m + 

l)T,x+(n±l)T) > q q (x). Recall r(t) = sup{x € R : u(t,a;) > 0}. If Nm > am, 
then r{mT) > amT, and furthermore, r(t) > (ara-l)f-L, (m— 1)T < t < mT. 
Hence from Lemma 12.71 we get a lower bound 

v e > v(e 2 ) - 10| log(l - p)| _1 = v(e 2 ) - 10| loge|~ 3 . (2.78) 

3 Comparison 

We now state precisely the comparison theorem we are using. Let U denote the 
set 

U = {(i, x) G [0, T] x R : x < vt} (3.1) 

for some v > 

Proposition 3.1. Assume |J.7| ), (CO|). 

i. Suppose that gift, u) < f{u) are Lipschitz functions, u G R, t > and 
initial data vq(x) < uq(x), x G R are given. There exists a probability space 
(QjJ-'jP) on which there are white noises W,W\, a solution u to a solu- 

tion v to 

[ d t v = dlv + g(v) + ea(v)Wi, (t,x)eU, ^ 
\v(t,x) = 0, (t,x)tU, [6 - Z) 

with v(0,x) = Vo(x), and satisfying 

u(t, x) > v(t, x), x G R, t > 0. (3.3) 

2. Fix possibly random -measurable uq G C. Suppose that g(u) is the Lipschitz 
function and there is also another initial data uq = vq G C cxp such that 

uo(-) < uo(-), a.s., 

and uo(x) = for all x > 0. Then there exists a probability space (£l,J-,P) 
on which white noises W, W, Wi,W2 and a vector of processes (u,u,v,w) are 
defined and satisfy the following properties: 

(i) u is a solution to 

(ii) u is a solution to M.l]) in C cxp starting at uq with W replacing W ; 
(Hi) v is a solution to \3. 2\) in C cxp starting at vq = uq where U is of the 

form iTO) : 

(iv) W<i is a space-time white noise independent ofW\; 



20 



(v) Up to time 

t = inf{t > : 3x, a^(t, x) =: a 2 (u) - a 2 {v) = 0, x) 7^ x)}, (3.4) 
u — v < w and u<u a.s., (3-5) 



where 
and 



d t w = d x w + w + ea w W 2 + S x - Vcom tA, (3.6) 



A(t) = - J (v(t,x)-v(0,x))dx+ J J g(v(s,x))dsdx+e J J a(v(s, x))W 1 (dsdx). 

Proof. 1. Assume first that er is globally Lipschitz. Then the proof goes essen- 
tially along the lines of Theorem 3.1 of |MP92| . One approximates the solutions 
by lattice versions as in (|1.15|) . for which the ordering is elementary. Then one 
shows the ordering is preserved in the limit. Because one has strong uniqueness 
it means the solution of the SPDE's are ordered in the desired way. Now sup- 
pose we do not have the strong uniqueness. We construct a sequence of Lipshitz 
er(") converging unifomly to a and consider the sequence of solutions 
corresponding to a^ n \ It is a standard to check that the sequence is tight and 
any weak limit point satisfies our equations. Since comparison is satisfied for 
each n it also holds in the limit. 

2. There exists a probability space with a noise W and a pair of independent 
noises W\ and Wi such that u solves (jl.ip . v solves (|3.2p and w solves 

d t w = d 2 x w + f(w + v)- f(v) + a w W 2 + 5 x - Kt A, (3.7) 

where crj = \o~ 2 (y + w) — a 2 (v)\, w is non-negative and 

u<v + w. (3.8) 

The construction of such a tripple (it, v, w) is fairly straightforward. One con- 
structs a sequence of approximations to (jl.ip . (|3.2p and (|3.7p for which the 
ordering correspondent to (|3 - 8[) is elementary. Then one takes a limit to get 
solutions to (jl.ip . (|3.2p and (|3 . 7|) and shows that ordering is preserved in the 
limit. By this way one gets that the the unique weak solution u to (jTTTJl is 
bounded from the above by v + w where v, w are some solutions to (|3.2p and 
(|3.7p respectively with independent white noises W\,W 2 - 
Define 

u = v + w, t < ti (3-9) 
8 t u = dlu + f{u)+ca(u)W,t>Tx, (3.10) 

where T\ is defined similarly to r: 

ti = inf{t > : Ete, &%,(t,x)=0, (v + w)(t, x) ^ v(t, x)}. 
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It is easy to see that u is a solution to the equation which u is supposed to 
solve, and hence we can set u = u and t\ — r. To show that u and v indeed 
belong to C cxp one can use for example the methods of proof of Theorem 1.2 
from [MPS06]. Now let us construct w satisfying (|3.6p such that 

w < w, out < T. 

Let w be a solution to 

d t w = d*w+(w + w)-(f(w + v)-f(v)), t<r. (3.11) 

As the drift term is non-negative we get that w is non-negative. Now define 
w = w + w and it is easy to check that it satisfies (|3.6|) and we are done. □ 



4 Large deviations 

We now present a fairly standard type of large deviation result which covers the 
estimates we need both in the upper and lower bounds. We need some notation. 
Let g(s,y,t,x) and r](x,y) be deterministic, and 

T b ,T = {(t,x):te [0,T],x-vt£ [6-1,6]}. (4.1) 

For (t,x) and (t',x') in T btT let 

d((t, x), {t',x')) = \x' -x\ + \t' - t\ x ' 2 . (4.2) 

Define 

„, ,s I f™[9(s,y,t',x') - g(s,y,t,x)] 2 rj(s,y)dsdy 
B{g,r),b)= sup — ^ . (4.3) 

«,*),(*, *>)er biT d((t,X), {t',x')) 

Lemma 4.1. Let g(s, y, t, x), r)(s, y), T b ^T and B(g, 77, b) be as above and a(t, x) 
nonanticipating with 

W(t, x)\ 2 < r)(t, x), (t, x) e [0, T] x R (4.4) 

almost surely, and define 

Z(t,x)= [ f g(s,y,t,x)a(s,y)W(dsdy). (4.5) 



There exist C r^g| , C r^~^| < oo such that if T > 1 and 

^Cq^TB 1 / 2 (4.6) 
Then, with = y/d(l + |log 2 d|), 

\Z{t',x')-Z{t,x)\ 



P sup 

\{t,x),{t',x')eT h , T ®{d({t 



(4.7) 
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Proof. Let Q n be the vertices of an afhne lattice with edges £ n parallel to the 
boundaries of Tb,T and with edge lengths 2~™ in the (1,0) direction and with 
vertical component 2~ 2n in the (l,v) direction. Let Q — U^L <?„ and £ — 

Given (t,x), (t',x') e Q with 

min(d((t, as), (*', x')), 1) € (2-"°, 2-<"°- 1 )] (4.8) 

there exists a path between them using edges from £, which uses only edges 
from £ n with n > uq, and uses at most T edges from any given £ n . 

For e = (p, q) £ £ n , write Z e := Z(p) — Z(q) and d e — d(p, q). By standard 
ltd calculus, 

£[exp{ 7 Z e }] < exp{i 7 2 d e 6}. (4.9) 
Let an = (10V2T)- 1 (n+ l) 1 /^-"/ 2 and 

A = {^e < a n i}. (4.10) 

By Chebyshev's inequality 

P{A%) < exp{i 7 2 4S - ja n £}. (4.11) 
Optimising the inequality over 7 gives 

P(^) < cxpl-Va 2 ^ 1 ^- 1 }. (4.12) 
Let A = n"=on e e5„A. On A, 

OO 

\Z(t',x')- Z(t,x)\ < TiJ^an (4.13) 

n=n 

< £(n + l) 1/2 2- no/2 (4.14) 

< l<$>(d((t,x),(t',x'))). (4.15) 

Here we use the fact that for n a > 0, E"=„ +i ™ 1/2 2~™ /2 < 10(n + l) 1/2 2-™°/ 2 . 
Now we have 

00 00 

p(A c ) < e E p (^) ^ E i^i^)- ( 4 - 16 ) 

n=0e££„ n=0 

It is simple to check that \£ n \ < 2T2 3 " and d e < 2~ n + 2 . From gUg) then, 

F(^ c ) < 2Te- 2 "° £2T " e "(l - e 3lo g 2-2-"/»r-»B- 1 j (417) 

which gives (|4. T[) as long as (|4.6[) holds. Since Z(s, y) is continuous, it is enough 
to check the bound on dyadics, and hence this completes the proof. □ 



23 



In order to apply Lemma we need a bound on (|4.3[) . This is provided by the 
next lemma. The lemma will only be applied with the A defined in (|2.14|) , but 
it is true for other A satisfying (|4.20p . 

Lemma 4.2. Let 

g(s,y,t,x) = e a ('- s )l(0 < s < t)G v (s,y,t,x) (4.18) 

and 

V(s, y)<(l + \v- vs|) exp{A|t/ - vs|}. (4.19) 
Then, there exists C ^.21\ ^ 00 suc ^ that for any b < (i) if T < 1 and a < 1, 
or (m) if T > 1 and 

A 2 

vA - — > 2a, (4.20) 

Z? /rom {4-.Sfy satisfies 

B<C^(l + \b\)e^. (4.21) 
TTie same a/so Zio/rfs m case (%) ?/G v is replaced by G v ,l- 

Proof. The only statement that is not elementary is (ii) . We have to estimate 



[g{s, y,t + h,x + z) - g(s, y, t, x)} (1 + \y - vs|) exp{A|y - vs\}dsdy 

<C{l + \b\)e m [z + h^ 2 } (4.22) 

for t £ [0, T], x - vt E [6-1,6], with g(s, y, t, x) in (|4.18[) . First of all, note that 
we can express G v in terms of Go, which in turn can be written explicitly in 
terms of the heat kernel G (see (U.13[) ) using reflection; 

G v (s, y, t, x) = e -i « a: - vt )-to- v;i ))-T(*—)G , (s, y - vs, t, x - vt), (4.23) 

G (s,y,t,x) = G(s,y,t,x) - G(s,y,t,-x), x,y<0. (4.24) 

After change of variables, the left hand side of (|4.22[) becomes, with x' = x — vt 

and 7 = z — vh, 

{lo^t+he-i^'-y+^-^-^Gois, y,t + h,x'+ 7 ) 

te -W-y)-^-s) G() ( Siy ^x'))\\ + \y\)e x ^dsdy. (4.25) 



LQ<s<tt 



with a = ^ — a. Estimating the two pieces of the right hand side of (|4.24| by 
using that the square of the sum is bounded by twice the sum of the squares we 
see that f|4. 25[) is bounded by the sum over l = ±1 of 

r , exp{-i(z' - y + 7 ) - a(s + h) ^ffl' } 

/ [ 1 -h<s<t , 

y <oJ v Vs + h 

VS / 47T 
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Note that we have also changed variables t — s s. Changing y i— > y + x' and 
rearranging a little this becomes {A\ + Ao\x'\)e~ where Ai is (47r) _1 times 
the sum over i = ±1 of 



dsdy 



//e— (1+lW(l ._ 1 ii:^p- l0S . s ,f=) 

(4.27) 

with 7' = exp{^ — a/i}. Consider the t = +1 term. We estimate 

, («-T) 2 _„ 2 (H-T) 2 _y 2 

/-|/g 4(s + h) g — \ 2 /g 4(s + h) g — \ 2 



-h<s<t , i-0<s<t 



) < 3 • lo< s <t 



ys + h V s ' v v s + h v >' 

,„ (y-7) 2 , , N „ (y— r) 2 

+3 • l_/,< s <o — 7 h 3 • lo< s <t — ; • (4.28) 

- - s + h ~ ~ s + h 

There are also three analogous terms corresponding to 1 = — 1. All six terms are 
estimated by explicit computation. Since it is very tedious, we present only the 
worst case which is the first term on the right hand side of (|4. 28[) with i = 1. 
Call P = v - A. 

Lemma 4.3. For h, 7 G (0, 1), a, <E (0, 3) and 

a = 2a - /? 2 > 0, (4.29) 
i/iere exists a < 00 such that for all t > 0, 

2 



" 2QS+ft (i + \v\) ( - < q^IM + H 1/2 ] 

(4.30) 

Proof. The left hand side is bounded by a constant multiple of Ii +I2 + Ia(7 2 — 
fc,) + h(h) + I 4 where 

Ii = f fe-^y(l + \y\) ( e -^- - e -jt) dyds 



!«! 



g 4( 


,+h) 


V 


S 




y 2 


e 4 < 


,+h) 




±h 




v 2 



2 



iv-ir 



I 3 (a) = / / e-^y{\ + \y\) I + I dyds 



0V( 7 2 -/i) 



- 2as+ , y (±tMl (e-^-e^^dyds 
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In the proof C will denote any finite constant, possibly depending on a,f3 and 
q. Its value will change from line to line. 

Estimation ofI± < Ch 1 / 2 . By the mean value theorem, there exists 9 G [0, 1] 
such that 

I 1 = —h 2 [ I ' e-^+^il + ^s^e-^myHs + ehy^dyds. (4.31) 
16 J h J 

Since we are integrating over s G [h, t] , we have s < s + 9h < 2s and therefore 
Ii < Ch 2 J* J e~ 2as+ P y e-£(y 4 + \y\ 5 ) S - 5 dyds 

= Ch 2 J e- qs s- 5 J e~ ^y- 2 ^ 2 (y 4 + \y\ 5 )dyds. 

Recall that \a + b\ n < 2"- 1 (|a|" + |6| n ). Thus y 4 = ((2s)iz + 2s(3) 4 < 2 5 z 4 s 2 + 
2 7 s 4 /3 4 , \y\ 5 = \{2s)^z + 2s(3\ b < 2^z b S i + 2 9 s 5 /? 5 . After the change of variables 
change variables to z = (y — 2s/3) / V2s and integrating we can bound the last 
term by 



/■t 

Cti 



f e- qs (s~i + s- 2 + s~i /3 4 + s? /3 5 ) ds (4.32) 

where C is a universal constant. Since s > h in the region of integration and 
h<l, this is bounded above by 



Ch 2 J e- qs (h~h- 2 + s- 2 + h-^13 4 + /i^/3 5 ) ds 



(4.33) 



The estimate then follows from jte qs s 2 ds < J, s 2 ds + f*e qs ds < \+q 1 . 

Estimation 0/I2 < Cbh . By the mean value theorem, there exists 9 e [0, 1] 
such that 

h = ~h 2 J J e - 2as+0v {l + \y\)e-^Kj(s + 9h)- 3 dyds 

< Ch 2 J* J e - 2 « s +^(l + \y\)e~ *&7 s - 3 dyds 
ft 

= Ch 2 e- qs s-i{^ +2s^ +2^s(3)ds 
Jh 

The estimate then follows from f h e~ qs s~'i (l+si+s)ds < ^s~^ + s~ 2 + s~i^j ds < 
Ch~i. 

Estimation of h(~f 2 -h)+I 3 (h) < C[hi+j}. This just uses j {l+\y\)e- {jL ^ dy < 

f(l + a + \y- a\)e~ i]L 2^dy < C((l + a)r? + r) and Aa - 1 > 0. 

Estimation of I4 < C7. First we change variables to r = s + h, with h 
constant, and use the fact that 2ar > and t < 00 to see that 

I 4 < e 2ah f f I(r, y)dydr (4.34) 

Av 7 2 J 
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where 

I(r,y) = r _1 (l + \y\) " e 2 e 4 — e 4 >- e 

< 2e -( 2 ^)V- 1 (l + | 2 ;|) e - iJt ^(e^-l) 2 (4.35) 
+2e _ (2a+ £ )v _ 1(1 + |y|) / _ e -i^^ 

=: 7 X (r, y) + h(r,y) 

We have (e 2 ^ — l) 2 < C7 2 and furthermore J(l + M)e _ ~ "2/ ' < /(l + 
7 + /3r + \y - 7 - / 3r|)e~ < "-' 3r) d y < C[r 1/2 (1 + 7 + fir) + r] so, 

I 1 {r,y)dydr < C7 2 . (4.36) 

Similarly we bound / ^(t", 2/)<^2/ by 

2 A e -( 2a+ ^ )rr -x (1+ ^ + 2 + |y _^_2 |)(e -i-^_ e -i^ )2 ^ (4 37) 
j z z 

We need two standard estimates: For r > 0, 

( e ~ ij T? 1 - _ e ~h) 2 dy < Cr'^-i 1 (4.38) 
|y - Jl(e-^ ~ e-^) 2 dy < C(r* 7 + 7 2 ) (4.39) 
Changing variables in (|4.37[) to z = y — fir and using (|4.38|) , (|4.39|) we have 

I 2 (r,y)dydr < C 7 (4.40) 



which gives the required estimate for I4 and completes the proof of the lemma. 

□ 

This gives us a lemma which controls the large deviations on a long time 
interval. 

Lemma 4.4. Let 

Z(t, x) = J J e- {t - s ^G v { S , y, t, x)a(s, y)W(dsdy) (4.41) 

where a is nonanticipating. Assume 

a 2 {t, x) < F(x - vt) + 3e~ x{x ' vt l (4.42) 
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Then for T be as in &2.50\) , 

P(3(t, x) G [0, T] x R, x < vt : Z(t, x) > e~ l e~ x( - x - yt) ) < 1/16 (4.43) 
and with M as in fa2.45\ ), 

P( sup \Z(t,x)\ > e - V/160) < 1/32 (4.44) 

0<t<T, vt-M<x<vt 

Proof. The left hand side of (|4.43[) is bounded by 

oo 

^ P(3(t, a:) G r_„, T : Z(t, x) > e'V") (4.45) 

n=0 

where T_ n x is defined in (|4.1|) . Applying Lemma |4. II for each n with 



£ = e- 1 e Xn /2 v / TlogT (4.46) 



using t' = and $ < 2V? 1 log T and ([4~2~Tj) . we obtain that ([4~45l is bounded 

by 

OO 

4| loge| 4 exp{-A(n + FtV"} (4.47) 

n=0 

with A = e _2 Cr^7j| loge|~ 12 | logloge _1 | _1 as long as 



I = e-V^VTbgT > C^TB 1 ' 2 (4.48) 

From (|2~T4) . A— l>lif<5> 0.2. Hence (n + l)- 1 e An > n + 1 for all n > and 
(|4.47[1 is bounded above by 

oo 

4|loge| 4 ^exp{-A(n + l)} < 8|loge| 4 e~ A (4.49) 

n=0 

as long as A > log 2 which bounded by 1/16 for e < eo- □ 

5 Proof of Lemma 12.21 

Recall that v and q are the solution of (|2.31[) and (|2.ip with initial data F as in 
UM- 7, L, T, M are as in (12^471) . dCT l. (I2~2^D . 

Lemma 12.21 is basically a result about how the stochastic perturbation of a 
partial differential equation (|2.31|) stays close to its deterministic version. Such 
theorems are fairly standard, but we need to stay close on a fairly long time 
interval [0, T] where T = 0{\ loge| 4 ) as required by ([2~46]) . First of all, let v be 
the solution of 

d t v = dlv + f(v) + ea(v)Wi, x<vt, , . 

v(t,x) = 0, x > vt, 1 ' 
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with v(0, x) = F(x) with v > v by the comparison theorem, and let g be the 
solution of (|2.2ip . It suffices to show that 

P (Bt € [0,T] : u(t,x) > x) + 3Ae~ A(a; - vt) for some x e R\ < 1/16. (5.2) 

We have 

d t (v - g) = d 2 x (v - g) + f(v) - f(g) + ea^W, x < vt (5.3) 
and v — g = on a; > vt and t = 0. One easily checks that 

f(v)-f(g)<2~(v-g). (5.4) 
Using the same ideas as in the proof of Proposition 13 . II we will show now that 

v — g < y (5-5) 

with y a solution of 

d t y = dly + 2-y + ea(v)W 1 . (5.6) 

on x < vt, and y = otherwise. To prove this define w to be the solution to 

d t w = d 2 x w + 2-(v-g)~ (/(«) - f(g)) - w. 

on x < vt, and w = otherwise. Note that by ([5"3)l . 2—(v—g)-(f(v)-f(Q)) > 
and hence w > 0. Now define 

y = u — £? + w, 

and by trivial calculations we get that y satisfies (|5.6p and since w > 0, (|5.5p 
follows. 

Using the integrating factor e* we obtain 

x) - g(t, x) < eZ(t, x) + 2 (5.7) 
where, with G v (s,y,t,x) as in (|2.32|) . 

Z(t,x)= J J e^ t - s) G v (s,y,t,x)a(s,y)W 1 (dsdy) (5.8) 

So it suffices to show that 

P(3(t,x) £ [0,T] x R,x < vt : > e -i e ->(*-v*)) < i/i 6 . (5.9) 

Note that when we do this we can assume without loss of generality that 

a 2 (t, x) < g(t, x) + ?>e- x{x - vt \ (5.10) 
For if v is a solution of (|3.2[) with cr 2 ({i(t,x)) replaced by 

5 2 {i,t,x) =min ( f r 2 (i;(t,x)),^,x)+3e- A ( a: - v *^ 

then 5 = w up to time 

f = inf{i > : v(t,x) > g(t,x) + 3e~ A(x - vt) for some x € R} 

and hence it suffices to prove (|5.9|) under (|5 . 10|) . The result now follows from 
Lemma 53] □ 
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6 Proof of Lemma 12.31 



i. of Lemma 16.11 is actually an upgrade of Lemma 12.21 which is a bit stronger 
than Lemma 12.31 ii. is similar to i. It is needed in Section [8] to control the 
maximum of w. Recall F from (|2.24[) and define 

( 0, x > 0, 

F(x) = l m*/40, xe(-M,0], (6.1) 
[ F[x) + 3Xe- Xx , x < —M, 

Lemma 6.1. i. Let v be the solution of \2.31\) with initial data \2.2J$ . Let 
q be the solution of \2.21\) with the same initial data. Let 7, L, T, M be as in 
(245^ - KM) . Then 

P (3t € [0, T] : v(t, x) > F(x - vt) for some x € R) < 1/8. (6.2) 

ii. Suppose that u satisfies < u(0, x) < F_(x) from 16.1]) and 

d t u < d 2 x u + u + ea(t, x)W (6.3) 

with 

o~(t, x) 2 < 3F(x ~vt — L). (6.4) 

Suppose M , 7 are as in {2.^5^ , f£.^7| ). Then there is a C rgy^ < 00 such that 

P( sup u(t,x)>u*/lO)<Cr K n 1 (6.5) 

0<t<3, -l<x-vt<l 

Let us make a few remarks before we start the proof of the lemma. Note that 
part (i) of the lemma is only stronger than Lemma 12.21 in the region x — vt £ 
(-M, 0] where g(t, x)+3Ae" A(;r " vt) ~ 3. We need this to be able to control o 2 (v) 
from below in the regions where v is small, that is in the region x~vt E {—M, 0]. 
If we will show that with high probability v is small in that region, then we will 
be able to use (| 1 . 5[) to control cr 2 (v) from below there. Another remark deals 
with coefficient 3 in (|6.4p . This coefficient appears in (|8 . 1 5f> after which we use 
Lemma 16.11 

Proof, i. Note also that by the same argument as that at (|5.10[) we can assume 
that 

a 2 (v(t,x)) <F(x-vt). (6.6) 
Let Af be the vertices of an affine lattice in 

T = {(i, x) : < t < T, vt - M < x < vt} 

with edge length ae between nearest neighbour vertices. From (|2.41[) we have 
that for (x, t) g T and therefore for p S Af, 

E[v(p)]< g<s 2 Me M . (6.7) 
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By Markov's inequality, 

P(v(p) > u*/40) < A0Me M e 2 /u*. (6.8) 

We can estimate 

P ( sup v(p) > 1*740] < P ( v (p) > u 7 40 ) ( 6 - 9 ) 



pert 



and since |r| < 2MTe- 2 Q- 2 , if a > 80M 1 / 2 e M / 2 (u*)- 1 / 2 then 



P sup v(p) > u*/40 < 1/32. (6.10) 
\peAf J 

Hence to prove the lemma it suffices to show that if a < C'jgTjgj 6 ^* then 



P sup \v(t',x')-v{t,x)\>u*/40 < 1/16. (6.11) 

I (t,x),(t',x')er J 

\ |x'-x| + |t'-t|<ae / 

Divide T into T intervals of length 1, T l = Tn{i < t < i + 1}. On {i < t < i + i}, 
v is the solution to 

v(t,x)= J G v (i,y,t,x)v(i,y)dy + J J G v (s,y,t,x)f(v(s,y))dsdy + eZ(t,x). 

(6.12) 

where 

Z(t,x) = J J G v {s,y,t,x)<j(v,s,y)W 1 {d S dy). (6.13) 
From (|6.6[) . assuming > 

K*',0 - w(t, x)\ < Qj. + n 2 + n 3 + e\Z(t',x')\ + e\Z(t, x)\. (6.14) 

where 

O x = 3 ( \G v (i iy ,t\x')-G v {i,y,t,x)\e^y-^dy, (6.15) 



^2 = 3 J J' \G v (s,y,t',x')-G v (s,y,t,x)\e- x( -y- vs Usdy, (6.16) 
Q 3 = 3 y J G v (s,y,t' \x')e- x{y ~ vs) dsdy. (6.17) 

So the result follows from (|4.44p and the elementary fact that there exists c < oo 
such that if \t - t'\ + \x - x'\ < cu* , A, v S [1, 2] then 

fii + fi 2 + ^3 < u*/80. (6.18) 
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ii. From 

„<*,*)<»./<;«,,„*, x) EW y + «•*«,.) (6.19) 



where 



Z(t,x)= j J e~ s G{s 1 y,t 1 x)a(s,y)W{dsdy) (6.20) 

From the definition of F_ it is clear that we can choose M so that for all < t < 3 
and — 1 < x — vt < 1, 

J e 4 G(0, y, t, x)F(y)dy < u* /20. (6.21) 

so the result follows from the following large deviation estimate whose proof is 
elementary as it only has to hold on time intervals of order 1: There exists a 
C|gj3j < oo such that for 7 is as in (|2.47|) . 

P( sup ee*\Z(t,x)\>u*/20)<Cir* r (6.22) 

0<t<3, -l<z-vt<l 

□ 



7 The critical mass 

The following elementary computation identifies the critical mass for survival. 
Lemma 7.1. Let W be a white noise and w(t,x) be a positive solution of 

d t w = d 2 x w + bw + *&\/wW , w(0,x)=w o (7.1) 
where is adapted with $ > $0 for some nonrandom -&Q > 0. Then 

P(w{t) = 0) > 1 - e h % 2 r 1 E[ j w (x)dx] (7.2) 



Proof. By considering w — e bt w we can assume without loss of generality that 
b = 0. If 

dt<t> = d 2 x 4> - tiltf (7.3) 

then 

exp{— J <j){t — s, x), w(s, x)dx)} (7-4) 

is a supermartingale in the s variable on [0, t]. The solution of (|7.3| with 
4>(0, x) = n is <p(t, x) — ($q< + n -1 ) -1 and hence 

.E[exp{-n / w(t, x)dx} \ T ] = exp{-(0§t + n" 1 )" 1 / w(0, x)dx}. (7.5) 

Taking n — > 00 we get 

P(/ w(t, x)dx = 0) = E[exp{-^ 2 t- 1 / w(0, x)da;}] (7.6) 

and the lemma follows from e~ x > 1 — x. □ 
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The next lemma is needed to control the support of such aw in short time 
intervals, in terms of the immigration. Note that we will only have to have 
reasonable control, and the actual scale are not critical here, as it is in the 
previous lemma. 

Lemma 7.2. Let W be a white noise and w be a solution of 

d t w = d 2 x w + bw + d^/ulW + dn, < t < 1 (7.7) 

with w(0,x) = and let tp be the hitting time of (— r, r) c ; 

tP = mf{t > : supp(w(t)) n (-r, rf ^ 0}. (7.8) 

Suppose that fi is a positive adapted measure on [0, l]xR with support in [0, 1] x 
(— r/2, r/2), and i9 is adapted and i9 > i?o > 0. Then, letting M = J J Q /i(dtdx), 

P(ip > 1) > 1 - 100r- 2 e b $Q 2 E[M}. (7.9) 

Proof. By considering w = d^e~ ht w we can assume without loss of generality 
that 6 = and i?o = 1- Furthermore, by symmetry it is enough to prove the 
result when ip is the hitting time of (— oo, — r] with a constant 50 instead of 100 
on the right hand side. First let us consider the case G J-q. Note that for 
any 5 > 

g s (x) = 12{x + r + S)- 2 (7.10) 
satisfies d 2 gs — <?f/2 for x > — r — 5/2. Let 

Xg(t) =exp j— J gs(x)w(t,x)dx + J f* gs{x)fx{dsdx)\ . (7-11) 
Then X$(t A ip) is a submartingale. In particular, 

E[X S {1 A 0) | T ] > E[X S (0) | T ] = 1. (7.12) 

Let us assume temporarily that /j,(t) G JF , < t < 1. Since g > 0, and 
3 < 48r~ 2 on supp(/i) we have M gs (l Aip) < M 9S (1) < 48r" 2 M, 

E[exp{-Jg s (x)w(lAip,x)dx} \ T a ] > exp {-48r- 2 M} . (7.13) 

Note that 

P(V> > 1) > liminf^o £ [exp {- J g 5 (a;)w(l A -0, £)ete} | .Fo] (7.14) 

which proves the lemma when fj,(t) € .Fo. 

For the general case, note first that (|7.7[) has the property that if W\ and W2 
are two solutions with measures fi\ and fj,2 and independent white noises W\ 
and VT2, then iwi + u>2 is a solution with measure /ii + /i2- 

We construct a probability space on which we have this setup with adapted 
1} and Hi(dx,dt) — fii(dx)St i (dt). Let ipi and ip2 be the corresponding hitting 
times of (— 00, — r]. From (|7.13[) and (|7.14p . conditioning on F ti instead of Tq 
we have P(tpi > T) > exp{-48r~ 2 Mi} for i — 1,2 with Mj = / fii(dx) and 
hence P(^i A ^2 > 1) > £ , [exp{-48r" 2 (Mi + M 2 )}]. 

A finite induction then gives the result for fi(dtdx) = ^2^=0 l x n{dx)5t TL (dt) 
with /i n (dx) G Tt n - We can then take limits to obtain the result for all adapted 
positive measures fx. □ 
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8 Proof of Lemma 12.41 



We will solve (13.6[) iteratively, on short time intervals of length 1 and show that 
we can kill the mass of w on each interval separately. The reason to do this is 
that the noise in l|3.6p , which is needed to kill the mass, is only of the correct 
order near x — vt where v is relatively small. So one has to show that the 
mass vanishes quickly, before the front moves ahead, and the noise is no longer 
available. We will do all our bounds on the event {v(t, x) < F_{t — vt)} where 
F_ is defined in (|6.1[) . To be more precise, define 



T„=mf{i>0 : v(t,x) > F(x - vt) for some ieR}. (8.1) 

Let 

v(t, •) = v(t A t v , •). 

Let Wi,k ■> k = 1, 2, . . . be a sequence of independent white noises which are 
also independent of W\. We construct a sequence of processes Wf., k = 1, 2, . . . 
by solving 



d t w k = d 2 x w k + ||/||upWfc + e<7kW 2 ,k + S x ^ vt A k k-l<t<k + l, 
w k (t,-) = t = k-l, 



(8.2) 



and setting w k (t, ■) = for t € (fc — 1, k+ l] c . Here 

u\ = \a 2 (v + w k + Wk-i) - u 2 (v + w k -i)\ V a*w k (8.3) 

and 

A k = Al k ^ 1<t < k (8.4) 

is the creation term acting only on the first half of each time interval. To start 
things going we use the convention that iB—i — wo = 0. Define stopping times 

T k ,x = inf{t e (k - l,k + 1] : supp{w) fe (t)} <f_ (v(k - 1) - l,vfc+ 1)}, 

r k ,2 = inf{* e (fc-2, fe + 1] :w k (t,x) + Wk-i{t,x)+v(t,x) >u*/l0 

for some x £ (v(fc — 2) — 1, vfc + 1)}, 

oo, if w k {k + 1, •) = 0, 

k + 1, otherwise. 



Tk,3 



with the convention that the infimum is infinite if the set is empty, and supp{u>} = 
{x : w(x) > 0} is the support of a non-negative function w. 

Let t = Tk,i be the smallest r k ,i < T, if there is one. Otherwise let r = T. 
Note that up to time r A (k + 1) we have 

a\ = a 2 (v + w k + tDfc-i) - u 2 {v + w k -i). 

Hence, we can find a probability space on which there are white noises W2 and 
{W2,i , i = l,2} such that the solution w of (|2.38[) can be represented as 

00 

w(t,x)l(t < t A t v ) = ^2w k (t,x)l(k - Kt < (k + 1) At At,). (8.5) 
fc=i 
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For each k = 1, . . . , T let 

U k = {r fc ,i > k + 1} n {r fe! 2 > k + 1} n {T fe ,3 > A: + 1} (8.6) 
Recall 7 from (|2.47|) . We claim that 

P(n{ =1 U k ) > l-c^T. (8.7) 
This implies Lemma |2~3I for on C\^ =1 Uk H {t^ > T}, 

oo 

w(*,x)=^tD fc (t,a;)l(fc-l<i<fc + l), t<T. (8.8) 
fe=i 

Hence the left hand side of (|2.53[) is bounded above by 

i - p (nLi^- n K > T}) < Co7 r + 1/8, 

since P(r„ < T) < 1/8 by Lemma \6.1l and we are done. 

The rest of the section will be devoted to verifying (|8.7p . Clearly it suffices 
to prove that for each k = 1, . . . , T and i = 1,2, 3, 

P(T fe ,i < k + 1) < c o7 . (8.9) 

By Lemma l7.1[ applied to idk on the interval (fc, fe+1] where there is no creation 
acting on u>k, 

P(r k , 3 < k + 1) < C^i^e^EjAjk) - A(k - 1)] < 2C D7 , (8.10) 

whereC jgjQ) , = e^H Li p(a*) _1 . In the last inequality we used (|2.44|) . By the 
same reasoning, but using Lemma lTT^l instead of Lemma FfTTl there is a C jg \i\ < 
oo such that 

P(T fc ,i<fc + l)<C m7 . (8.11) 
It remains to prove (|8.9p for rfe,2- The rest of the proof is devoted to this. Define 

u k (t)=w k (t,x)+w k -i(t,x) + v(t,x), t€[k-2,k + l]. (8.12) 

Given that Wkik — 1, •) = Wk-2(k — 1, •) = there is a white noise such 
that 

<9 t u fe = dlu k + (w k + Wk-i) + f(v) + £<J(k)W( k ) , (8.13) 
where <r? fc N = o\_ x + a\ + u{v) 2 . Since f(v) < v, 

d t u k < d 2 x u k + u k + ea {k} W ik ) ■ (8.14) 

The inequality is meant as holding for the corresponding integral equation. Now 
we claim that on {r^i > k + 1} n {rfe_i,i > k}, we have, for t < Tk,2 

ff( fc )(*,ar) < (||a 2 || Lip + 2)F(a ; -vt-L). (8.15) 
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For we know that there a 2 (v(t,x)) < v(t,x) < F(x — vt) < F(x — vt — 

L). AlsO_al(t,x) < \\(J 2 \\upW k (t,x) < ||CT 2 ||Lip^l{(fc-l,fc+l]x(v(fc-l)-l,vfc+l)} < 

\\a 2 \\ Up F(x-vt~L). 
Note that 

P{rk,2 < k+l,n,i > k+l,Tk-i,i > k) < P[ sup u k (t,x) > u*/X0) . 

v(fc-2)-l<^<vfc + l) 

(8.16) 

We want to estimate this with Lemma f6.ll but we need (|6.4[) . which does not 
necessarily hold. But by the argument around (|5.10|) . in proving (|6.4[) . we can 
assume without loss of generality that (|8 . 15[) holds. In fact this is the place 
where it is clear the appearance of coefficient 3 in (|6.4| . Hence we can apply 
Lemma 16. II to obtain 

P(r k ,2 < fc+l,7*,i > k + l,T fc _i ; i > k) < Cg3j 7 . (8.17) 

Finally, 

P(T k ,2<k+l) < P(T k ,2<k + l,T kA > k+l,T k - h l>k) + P(T k -l,l<k) 

+ P(r M < fc + 1) 
< (Q631 + 2 QTOh, 

which completes the proof of Lemma 12.41 

9 Proof of Lemma 12.61 

We need a preliminary result of how a stochastic perturbation of a partial dif- 
ferential equation stays close to its deterministic version. Here the interval is of 
order 1, but the estimate needs to be precise. 

Lemma 9.1. Suppose that u and g are solutions of 

d t u = d 2 x u + f(u)+a(u)W (9.1) 

and 

d t Q = d 2 x Q + f(g) (9.2) 

on \x\ < L + vt, < t < T with u(t,x) — g{t,x) — on \x\ > L + vt, and 
ii(0,x) = p(0,x). Suppose that f is Lipschitz with constant K. Then 

\u{t,x) - g{t,x)\ < \Z{t,x)\ + \Z{t,x)\. (9.3) 

where 

Z(t, x) = K J J e K ( 4 - s >G v , L ( S , y, t, x)\Z(s, y)\dsdy, (9.4) 

Z(t,x) = J J G v>L (s,y,t,x)a(u(s,y))W(dsdy). (9.5) 
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Proof. Let D = u — g — Z . Note that D satisfies 

d t D = dlD + f(D + Q + Z)-f( Q ) (9.6) 

on \x\ < L + vt with D(t, x) = on |ac| > L + vt and t = 0. Now \ f(D + g + 
v) - f(g)\ < K\D\ + K\v\. Let D + and D_ be the solutions of 

dtD±=8£D±±K[\D±\ + \Z\], \x\<L + vt (9.7) 

with D± (t, x) = for \x\ > L + vt and D±(0, x) = 0. Note that D+ is a super- 
solution and D_ is a subsolution of (|9.6p . so that D_ < D < D + . Furthermore, 
Z3_ < < D + , so (|9.7I ) can be solved explicitly in terms of Z. We get 

D±(t,x)=±K f I e K ^G v , L ( S ,y,t,x)\Z( S ,y)\dsdy. (9.8) 



as desired. □ 
Now we continue with the proof of Lemma 12.61 Define 

g(t,x) = g(t,x) +rey/gg(t,x) (9.9) 

%) = (l + H)exp{(l-^)|x|} (9.10) 



Proof of Lemma \2.(k Let 



A=\uj: \u{t,x) - g(t,x)\ < er^Gg{t,x), < t < l,x £ R| . (9.11) 
In particular, on A, we have 

u(t,x) < g{t,x), 0<i<l,ieR (9.12) 



If we were to let a(t, x) — y // mm{o~ 2 (t, x), g(t, x)}, and u be the solution of 
(|2.66p with a replaced by <7, and A the analogue of A with u replaced by u, 
then P(A) — P(A). Hence in estimating P(A) we can assume without loss of 
generality that 

a 2 (t,x) < g(t,x). (9.13) 



/ is Lipschitz with constant 1, so by Lemma [9TTI 

\u(t,x)-g(t,x)\ < e\Z(t,x)\ + e\Z(t,x)\, (9.14) 

where Z(t, x) and Z(t, x) are is as in ()9.4j) and (19. 5ft with u replaced by u and 
^=1. 

Now note that there exists C jQ < oo such that for < t < 1, 

Q{t,x)<C^^e 2 1 h{x-vt + L). (9.15) 
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By Lemmas 14.11 and 14.21 with T = 1, and a = 1, and by (|9.15[) . we have, for 
some C |g 161 , < oo, 

P( rap \Z(t,x)\ > re^y/h^}) < 4exp{-C | ^ ra r 2 }. (9.16) 

a— vte[b~— 1,6] 

Furthermore there exists C jg Tjg < oo such that for < t < 1, 

e 2 7 /t(x - vf + L) < C ^jT\ Gg{t 7 x) x < vt + L (9.17) 
so summing b from to L and using (|9.17p . 

P(e|Z(l,s)| < reJgg(l,x) for all x e [v, v + L)) > 1 - 4Lvip{-Cgj^r 2 }. 

(9.18) 

Finally, it is not hard to check that there is a C jq jg[ < oo such that for 
< t < 1, 

t 

e t - s G v ,L{s,y,t,x)^/h&TIjdsdy < C[gj§[ ^h(x - vt + L), (9.19) 
and therefore we also have for some C jg 2Q[ < oo, 

P(e|Z(l,a)| <reJgg(t,x) for all a; G [v, v + L)) > 1 - 4L cxp{-C|^r 2 }, 

(9.20) 

which completes the proof of Lemma 12.61 □ 
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